IFIC/ 11-70 



Supersymmetric mass spectra and the seesaw type-I scale 

C. Arbelaez,* M. Hirsch,^^ and L. Reichert"'- 
AHEP Group, Instituto de Fisica Corpuscular - C.S.I.C./Universitat de Valencia 
Edificio de Institutos de Paterna, Apartado 22085, E-46071 Valencia, Spain 

Abstract 

We calculate supersymmetric mass spectra with cMSSM boundary conditions and a type-I seesaw 
mechanism added to explain current neutrino data. Using published, estimated errors on SUSY 
mass observables for a combined LHC-I-ILC analysis, we perform a theoretical analysis to identify 
parameter regions where pure cMSSM and cMSSM plus seesaw type-I might be distinguishable 
with LHC-I-ILC data. The most important observables are determined to be the (left) smuon and 
selectron masses and the splitting between them, respectively. Splitting in the (left) smuon and 
selectrons is tiny in most of cMSSM parameter space, but can be quite sizeable for large values 
of the seesaw scale, mss- Thus, for very roughly mss > 10^^ GeV hints for type-I seesaw might 
appear in SUSY mass measurements. Since our numerical results depend sensitively on forecasted 
error bars, we discuss in some detail the accuracies, which need to be achieved, before a realistic 
analysis searching for signs of type-I seesaw in SUSY spectra can be carried out. 
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I. INTRODUCTION 



The seesaw mechanism [1-5] provides a rationale for the observed smallness of neutrino 
masses [6-11]. However, due to the large mass scales involved, no direct experimental test 
of "the seesaw" will ever be possible. Extending the standard model (SM) only by a seesaw 
mechanism does not even allow for indirect tests, since all possible new observables are 
suppressed by (some power of) the small neutrino masses. ^ 

The situation looks less bleak in the supersymmetric version of the seesaw. This is es- 
sentially so, because soft SUSY breaking parameters are susceptible to all particles and 
couplings which appear in the renormalization group equation (RGE) running. Thus, as- 
suming some simplified boundary conditions at an high energy scale, the SUSY softs at the 
electro-weak scale contain indirect information about all particles and intermediate scales. 
Perhaps the best known application of this idea is the example of lepton flavour violation 
(LFV) in seesaw type-I with cMSSM ^ boundary conditions, discussed already in [14]. A 
plethora of papers on LFV, both for low-energy and for accelerator experiments, have been 
published since then (for an incomplete list see, for example, [15-27]), most of them concen- 
trating on seesaw type-I. 

Seesaw type-I is defined as the exchange of fermionic singlets. At tree-level there is also 
the possibility to exchange (Y=2) scalar triplets [4, 5], seesaw type-II, or exchange (Y=0) 
fermionic triplets, the so-called seesaw type-Ill [28, 29]. Common to all three seesaws is that 
for m,y ~ a/ Am\ ~ 0.05 eV, where Am\ is the atmospheric neutrino mass splitting, and 
couplings of order 0(1) the scale of the seesaw is estimated to be very roughly mss ~ 10^^ 
GeV. Much less work on SUSY seesaw type-II and type-Ill has been done than for type-I. 
For studies of LFV in SUSY seesaw type-II, see for example [30, 31], for type-Ill [32, 33]. 

Apart from the appearance of LFV, adding a seesaw to the SM particle content also leads 
to changes in the absolute values of SUSY masses with respect to cMSSM expectations, at 
least in principle. Type-II and type-Ill seesaw add superfields, which are charged under 
the SM group. Thus, the running of the gauge couplings is affected, leading to potentially 
large changes in SUSY spectra at the EW scale. In [34] it was pointed out, that for type-II 
and type-Ill seesaw certain combinations of soft SUSY breaking parameters are at 1-loop 
order nearly constant over large parts of cMSSM parameters space, but show a logarithmic 
dependence on niss- ^ This was studied in more detail, including 2-loop effects in the RGEs, 
for type-II in [31] and for type-Ill in [32]. Using forecasted errors on SUSY masses, obtained 
from full simulations [36, 37], the work [38] calculated the error with which the seesaw 



^ "Low-energy" versions of the seesaw, such as inverse seesaw [12] or linear seesaw [13], might aUow for 
larger indirect effects. In this paper we will focus exclusively on the "classical" seesaw with a high (B-L) 
breaking scale. 

^ "constrained" Minimal Supersymmetric extension of the Standard Model, also sometimes called mSugra 
in the literature. 

^ These so-called invariants can be useful also in more complicated models in which an inverse seesaw is 
embedded into an extended gauge group [35] . 
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(type-II and -III) scale might be determined from LHC and future ILC [39] measurements. 
Interestingly, [38] concluded that, assuming cMSSM boundary conditions, ILC accuracies on 
SUSY masses should be sufficient to find at least some hints for a type-II/type-III seesaw, 
for practically all relevant values of the seesaw scale. 

Seesaw type-I, on the other hand, adds only singlets. Changes in SUSY spectra are 
expected to be much smaller and, therefore, much harder to detect. Certainly because of 
this simple reasoning much fewer papers have studied this facet of the type-I SUSY seesaw 
so far. Running slepton masses with a type-I seesaw have been discussed qualitatively in 
[23, 24, 40, 41]. In [42] it was discussed that in cMSSM extended by a type-I seesaw, splitting 
in the slepton sector can be considerably larger than in the pure cMSSM. This is interesting, 
since very small mass splittings in the smuon/selectron sector might be measurable at the 
LHC, if sleptons are on-shell in the decay chain X2 ~^ ^^^^ ~^ ^^^^Xi [43]. 

In this paper, we calculate SUSY spectra with cMSSM boundary conditions and a see- 
saw type-I. We add three generations of right-handed neutrinos and take special care that 
observed neutrino masses and mixing angles are always correctly fitted. We then follow the 
procedure of [38]. Using predicted error bars on SUSY mass measurements for a combined 
LHC-I-ILC analysis, we construct fake "experimental" observables and use a x^-analysis to 
estimate errors on the parameters of our model, most notably the seesaw scale. We identify 
regions in parameter space, where hints for a type-I seesaw might show up at the ILC/LHC 
and discuss quantitatively the accuracy which need to be achieved, before a realistic analysis 
searching for signs of type-I seesaw in SUSY spectra can be carried out. 

The rest of this paper is organized as follows. In the next section we define the super- 
symetric seesaw type-I model, fix the notation and define the cMSSM. In section III we 
present our results. After a short discussion of the procedures and observables in section 
III A, we show a simplified analysis, which allows to identify the most important observables 
and discuss their relevant errors in section IIIB. Section IIIC then shows our full numerical 
results. We then close with a short summary and discussion in section IV. 

II. SETUP 

A. Supersymmetric seesaw type-I 

In the case of seesaw type-I one postulates very heavy right-handed neutrinos with the 
following superpotential below the GUT scale, Mg'- 

Wi = Wmssm + W,. (1) 

Here Wmssm is the usual MSSM part and 

= Ntyfij ■ Hu + lmMR,uN[ ■ (2) 

We have written eq. (1) in the basis where and the charged lepton Yukawas are diagonal. 
In the seesaw one can always choose this basis without loss of generality. For the neutrino 
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mass matrix, upon integrating out the heavy Majorana fields, one obtains the well-known 
seesaw formula 

= -^r-.^M-iy^ (3) 

valid up to order 0{m£,/Mji), m^) = -^y^- Being complex symmetric, the light Majorana 
neutrino mass matrix in eq. (3), is diagonalized by a unitary 3x3 matrix U [4] 

= U'^ ■ ■ U . (4) 

Inverting the seesaw equation, eq. (3), allows to express as [44] 



y = V2—\/Mr-R-^/^,-U\ (5) 

where the rhj, and are diagonal matrices containing the corresponding eigenvalues. R is 
in general a complex orthogonal matrix. Note that, in the special case R = 1, contains 
only "diagonal" products \fM^mi. For V we will use the standard form 



U = -S12C23 - Ci2S23Si3e*^ C12C23 - Si2S23Sl3e*^ S23C13 

\ S12S23 - Ci2C23Sl3e*'^ -C12S23 - Si2C23Si3e*'^ C23C13 / 



X 



/ e»"i/2 \ 
e^"2/2 

V 



(6) 



with Cij = cos 6ij and Sij = sin 6ij. The angles 612, O13 and 6*23 are the solar neutrino angle, 
the reactor angle and the atmospheric neutrino mixing angle, respectively. 6 is the Dirac 
phase and a, are Majorana phases. Since U can be determined experimentally only up to 
an irrelevant overall phase, one can find different parameterizations of the Majorana phases 
in the literature. 

Eq. (3) contains 9 a priori unknown parameters, eq. (5) contains 18. The additional 9 
unknowns encode the information about the high scale parameters, the three eigenvalues of 
Mr and the 3 moduli and 3 phases of R. 

B. cMSSM, type-I seesaw and RGEs 

The cMSSM is defined at the GUT-scale by: a common gaugino mass M1/2, a common 
scalar mass mo and the trilinear coupling Aq, which gets multiplied by the corresponding 
Yukawa couplings to obtain the trilinear couplings in the soft SUSY breaking Lagrangian. 
In addition, at the electro- weak scale, tan/3 = v^/vct is fixed. Here, as usual, Vd and are 
the vacuum expectation values (vevs) of the neutral component of Hd and H^i respectively. 
Finally, the sign of the /i parameter has to be chosen. 

Two-loop RGEs for general supersymmetric models have been given in [45]. ^ In our 
numerical calculations we use SPheno3.1.5 [47, 48], which solves the RGEs at 2- loop, in- 
cluding right-handed neutrinos. It is, however, useful for a qualitative understanding, to 



^ The only case not covered in [45] is models with more than one U{1) gauge group. This case has been 
discussed recently in [46]. 
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consider first the simple solutions to the RGE for the slepton mass parameters found in the 
leading log approximation [16, 22], given by 

(AM|),, = -^,{'iml + Al){Y^^^LY'^),, (7) 
{AMl)ij = 0, 

where only the parts proportional to the neutrino Yukawa couplings have been written. The 
factor L is defined as 

Eq. (7) shows that, within the type-I seesaw mechanism, the right slepton parameters do not 
run in the leading-log approximation. Thus, LFV is restricted to the sector of left-sleptons 
in practice, apart from left-right mixing effects which could show up in the scalar tau sector. 
Also note that for the trilinear parameters running is suppressed by charged lepton masses. 

It is important that the slepton mass-squareds involve a different combination of neutrino 
Yukawas and right-handed neutrino masses than the left-handed neutrino masses of eq. (3). 
In fact, since [Y'^'^^LY'^) is a hermitian matrix, it obviously contains only nine free parameters 
[17], the same number of unknowns as on the right-hand side of eq. (5), given that in principle 
all 3 light neutrino masses, 3 mixing angles and 3 CP phases are potentially measurable. 

Apart from the slepton mass matrices, Y'^ also enters the RGEs for at 1-loop level. 
However, we have found that the masses of the Higgs bosons are not very sensitive to the 
values of Y'^, see also next section. We thus do not give approximate expressions for m^Hu- 
For all other soft SUSY parameters, Y^ enters only at the 2-loop level. Thus, the largest 
effects of the SUSY type-I seesaw are expected to be found in the left slepton sector. 



III. NUMERICAL RESULTS 
A. Preliminaries 

We use SPhenoS .1.5 [47, 48] to calculate all SUSY spectra and fit the neutrino data. 
Unless noted otherwise the fit to neutrino data is done for strict normal hierarchy (i.e. 
f^ui = 0), best-fit values for the atmospheric and solar mass squared splitting [10] and 
tri-bimaximal mixing angles [49]. To reduce the number of free parameters in our fits, we 
assume right-handed neutrinos to be degenerate and R to be the identity. The seesaw scale, 
called mss below, is equal to the degenerate right-handed neutrino masses. We will comment 
on expected changes of our results, when any of these assumptions is dropped in the next 
subsections. Especially, recently there have been some indications for a non-zero reactor 
angle, both from the long-baseline experiment T2K [50] as well as from the first data in 
Double CHOOZ [51]. We will therefore comment also on non-zero values of ^13 = 6r. 
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SPheno solves the RGEs at 2-loop level and calculates the SUSY masses at 1-loop order, 
except for the Higgs mass, where the most important 2-loop corrections have been imple- 
mented too. Theoretical errors in the calculation of the SUSY spectrum are thus expected 
to be much smaller than experimental errors at the LHC However, since for the ILC one 
expects much smaller error bars, theory errors will become important at some point. We 
comment on theory errors in the discussion section. 

Observables and their theoretically forecasted errors are taken from the tables (5.13) 
and (5.14) of [36] and from [37]. For the LHC we take into account the "edge variables": 
{miiy^^^, (mig)lJ^ (m/,)^'J;, (m^ijcdge and (m^^Jthrcsh from the decay chain qi and 
X2 ^ ^^Xi [52-54]. In addition, we consider {mubjthresh, ("^r+r-) (from decays involving 
the lighter stau) and the mass differences A-^^ = mg — mi^, with z = 1, 2, A^^^o = m,^ ""^x? 
and A;-^^o = mj^ — m^o. Since m^^ — ~ m^^ ~ m^^ applies for a large range of the 
parameter space LHC measurements will not be able to distinguish between the first two 
generation squarks. The combined errors for an LHC+ILC analysis, tables (5.14) of [36], 
are dominated by the ILC for all non-coloured sparticles, except the stau. For us it is 
essential that both, left and right sleptons are within reach of the ILC. Also the two lightest 
neutralinos and the lighter chargino measured at ILC are important. The errors in [36] were 
calculated for relatively light SUSY spectra, thus we extrapolate them to our study points, 
see below, assuming constant relative errors on mass measurements. We will comment in 
some detail on the importance of this assumption below. Finally, we use the splitting in the 
selectron/smuon sector [43] as an observable: 

A(m--) = ~ (9) 

Here, mp^"" = |(me- + m^). The LHC can, in principle, measure this splitting from the edge 
variables for both, left and right sleptons, if the corresponding scalars are on-shell. In cMSSM 
type-I seesaw only the left sector has a significant splitting, we therefore suppress the index 
"L" for brevity. For this splitting [43] quote a "one sigma observability" of A(me^) ~ 2.8 
%o for SPSla. ^ For comparison, the errors on the left selectron and smuon mass at the ILC 
for this point are quoted as A{mg) ^ 1 %o and A(m^) ^ 2.5 %o, respectively [36]. 

The negative searches for SUSY by CMS [55] and ATLAS [56] define an excluded range 
in cMSSM parameter space, ruling out the lightest SPS study points, such as SPSla' [37] 
or SPS3 [57]. For our numerical study we define a set of five points, all of which are chosen 
to lie outside the LHC excluded region, but have the lightest non-coloured SUSY particles 



SPSla has only the edge in the right-slepton sector on-shell, see discussion fig. (3). 
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within reach of a 1 TeV hnear colhder. The points are defined as follows: 



Pi ^ (mo = 120, Ml/2 = 600,^0 = 0,tan/3 = 10) 

P2 ^ (mo = 120, Ml/2 = 600, Aq = 300, tan/3 = 10) 

P3 ^ (mo = 120, Ml/2 = 600, v4o = -300, tan /3 = 10) 

P4 -> (mo = 180, Ml/2 = 550, Aq = 0, tan/3 = 10) 

P5 (mo = 180, Ml/2 = 550,^0 = 300, tan /3 = 10) 



(10) 



All points have sgn{fi) > 0, masses are in units of GeV. Points P1-P3 lie very close to the 
stau-coannihilation line. We have checked by an explicit calculation with MicrOmegas [58- 
61] that the relic density of the neutralino agrees with the current best fit value of QcDAih"^ 
within the quoted error bars [11] for Pi. P4 and P5 have been chosen such that deviations 
from the pure cMSSM case are larger than in P1-P3, see eq.(7), i.e. to maximize the impact 
of the seesaw type-I on the spectra, see below. 

B. Observables and seesaw scale 

In this subsection we will first keep all parameters at some fixed values, varying only the 
seesaw scale. These calculations are certainly simple-minded, but also very fast compared to 
the full Monte Carlo parameter scans, discussed later. However, as will be shown in the in 
the next subsection, there is nearly no correlation between different input parameters. Thus, 
the simple calculation discussed here already gives a quite accurate description of the results 
of the more complicated minimization procedures of the "full" calculation. Especially, this 
calculation allows us to identify the most important observables and discuss their maximally 
acceptable errors for our analysis. 

In fig. (1) we show 



where A(mj) is the expected relative experimental error for the mass of sparticle i at the ILC, 
as a function of 77255. We remind the reader that we assume that A(mj) can be extrapolated 
to our study points. To the left results for Pi and to the right for P5. ml^'^^^^'^ is the value 
of the mass calculated in the cMSSM limit and m["^® the corresponding mass for a seesaw 
scale of mss- These latter values have always been calculated fitting the Yukawa matrix of 
the neutrinos at mss-, such that the best fit values of solar and atmospheric neutrino mass 
differences are obtained and niy^ = is maintained. As expected the departures from the 
cMSSM values then increase with increasing seesaw scale. Note that the lines stop at values 
of mss ~ (2 — 3) X 10^^ GeV, since for larger values neutrino Yukawas, which are required 
to fit the neutrino data, are non-perturbative. 

Significant departures with respect to the cMSSM values are found (with decreasing 
importance) for the following observables: left smuon mass, left selectron mass, mass of Xi, 
m^o and xt- We have checked that all other observables have much milder dependences 




(11) 
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FIG. 1: Calculated deviations of masses from their nominal cMSSM values as function of mss for 
the most important masses. To the left Pi, to the right P5. 



on mssj as expected. The smuon mass is more important than the selectron mass, despite 
the latter having a smaller predicted error, due to our choice of degnerate right-handed 
neutrinos in the fits. With this assumption the running of the smuon mass has contributions 
from Yukawas responsible for both, atmospheric and solar scale, while the selectron has 
contributions from the Yukawas of the solar scale only. The change in Xi and xt masses 
are small in absolute scale, but it is expected that ILC will measure these masses with very 
high accuracy. Also m/jo shows some mild dependence on rriss, but on a scale of an expected 
experimental error of 50 MeV [37], i.e. much smaller than our current theoretical error, see 
below. 

As the figure shows deviations from cMSSM expectations of the order of several standard 
deviations are reached for left smuon and selectron for values of mss above 10^^ GeV. 
Comparing the results for Pi (left) with those for P5 (right) it is confirmed that P5 shows 
much larger deviations from cMSSM. We have checked that results for the other points P2-P4 
fall in between the extremes of Pi and P5. Lines for P2 and P3 are nearly indistinguishable 
in such a plot, apart from some minor difference in the Higgs mass. 

In fig. (2) we show the calculated x^ as a function of mss for 4 different cMSSM points. 
Here, is calculated with respect to cMSSM expectations. To the left we show Xt including 
all observables, to the right Xt without the mass splitting in the (left) smuon-selectron 
sector. The figure demonstrates again that Pi (P5) has the smallest (largest) departures 
from cMSSM expectations. A non-zero value of Aq can lead to significant departures from 
cMSSM expectations. Determination of Aq from measurements involving 3rd generation 
sfermions and the lightest Higgs mass will therefore be important in fixing mss- 

Fig. (2) also demonstrates that A(me^) at its nominal error gives a significant contri- 
bution to the total x^- Thus, LHC measurements only might already give some hints for a 
type-I seesaw [42]. However, with the rather large error bars of mass measurements at the 
LHC it will not be possible to fix the cMSSM parameters with sufficient accuracy to get a 
reliable error on the value of mss- Unfortunately, also the accuracy with which A(me/i) can 
be measured at the LHC is quite uncertain. According to [43] such a splitting could be found 
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mss [GeV] m.ss [GeV] 



FIG. 2: Calculated ^ function of 77155 for 4 different cMSSM points. To the left: Total 
including all observables, to the right total XT without the mass splitting in the (left) 

smuon-selectron sector. Values quoted in the plots correspond to (7770, -^^1/21 ^o)- I11 sll points 
shown we choose tan /3 = 10 and ^ > 0. 




J013 IQl-l IQl'' 1013 Ifll-l IQlS 

mss [GeV] mss [GeV] 



FIG. 3: Calculated x^ for the observable A(777e^) as function of 77155 for different values of its error. 
To the left: Pi; to the right P5. 

for values as low as (few) 10~^ or as large as (several) percent, depending on the kinematical 
configuration realized in nature. Moreover, our points P1-P5 have heavier spectra than the 
ones studied in [43], so larger statistical errors are to be expected. 

Fig. (3) shows the relative deviation of A(me/i) for Pi (left) and P5 (right) for different 
assumed values of the error in this observable, relative to cMSSM. Here, a = 1, 2, 3, 4 means 
that we have multiplied the "error" quoted in [43] by factors 1,2,3,4. The deviation drops 
below one sigma for any value of mss shown for Pi (P5) when this error is larger than twice 
(six times) the nominal error. This implies that no hints for seesaw type-I can be found in 
LHC data if the error on A{mgfi) is larger than 5 %o (1.6 %) in case of Pi (P5). 

We should also mention that the actual value of A{mep,) is not only a function of mss 
and the cMSSM parameters, but also depends on the type of fit used to explain neutrino 
data. We have used degenerate right-handed neutrinos and m^^ = in the plots shown 
above. Much smaller splittings are found for (a) nearly-degenerate light neutrinos, i.e. 
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FIG. 4: To the left Br(/i — )• cy) and to right calculated function of mss for two different 

values of the reactor angle 

rriiy^ > 0.05 eV; or (b) very hierarchical right-handed neutrinos. We have checked by an 
expUcit calculation that, for example, for P5 and m,^^ = 0, Ax^ > 5.89 ^ for values of 77155 
larger than mss — 1-6 x 10^^ GeV from A(me/i) alone, whereas the same Ax^ is reached for 
m^^ = 0.05 eV only for 77155 ^ 7 x 10^^ GeV. Consequently, even though one expects that 
a finite mass difference between left smuon and selectron is found in cMSSM type-I seesaw, 
this is by no means guaranteed. 

Similar comments apply to the errors for the selectron and smuon mass at the ILC. For 
Pi (-P5) the departure of the left selectron mass from the cMSSM expectations is smaller 
than 1 a even for 77755 ~ 3 x 10^^ if the error on this mass is larger than 1.5%o (1%). For 
the left smuon the corresponding numbers are for Pi and P5 approximately 1.5% and 5%, 
respectively. 

Naively one expects LFV violation to be large, whenever the neutrino Yukawa couplings 
are large, i.e. for large values of 77755. That is, the regions testable by SUSY mass mea- 
surements could already be excluded by upper bounds on LFV, especially the recent upper 
bound on /i — 67 by MEG [62]. That this conjecture is incorrect is demonstrated by the 
example shown in fig. (4). In this figure we show the calculated Br(/i — )■ 67) to the left and 
the calculated (total and only A{mgfi)) to the right for 5 = vr and two different values of 
the reactor angle, 6*13 for the point Pi. For 6*13 = all values of 77755 above approximately 
77755 ~ lO^'' GeV are excluded by the upper bound Br(/i — 67) < 2.4 x 10"^^ [62]. For 
^13 = 6° nearly all values of 77755 become allowed. At the same time, this "small" change in 
the Yukawas has practically no visible effect on the calculated from mass measurements 
as the plot on the right shows. This demonstrates that SUSY mass measurements and LFV 
probe different portions of seesaw type-I parameter space, contrary to what is sometimes 
claimed in the literature. That one can fit LFV and SUSY masses independently even for 
such a simple model as type-I seesaw is already obvious from eq. (7): Even after fixing all 
low energy neutrino observables we still have nine unknown parameters to choose from to 



^ Ax^ > 5.89 corresponds to 1 cr c.l. for 5 free parameters. 
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178 179 180 181 182 100 200 300 400 500 

mo [GeV] Ao [GeV] 

FIG. 5: Calculated allowed parameter space for uiq, M1/2, tan/?, Aq and mss for 7 free parameters, 
P5 and mss = 5 x 10^^ GeV. For discussion see text. 

fit any entry of the left slepton masses independently. 

Fig. (4) also shows that non-zero values of 6*13, as preferred by the most recent experi- 
mental data [50, 51], should have very little effect on our parameter scans. In our numerical 
scans, discussed next, we therefore keep ^13 = unless mentioned otherwise. We will, how- 
ever, also briefly comment on changes of our results, when 6*13 is allowed to float within its 
current error. 

C. Numerical scans 

For the determination of errors on the cMSSM parameters and mss we have used two 
independent programmes, one based on MINUIT while the other uses a simple MonteCarlo 
procedure to scan over the free parameters. For a more detailed discussion see [38]. Plots 
shown below are obtained by the MonteCarlo procedure, but we have checked that results 
from MINUIT and our simplistic approach described above give very similar estimates for 
the x^ with MINUIT only slightly improving the quality of the fit. In this section we always 
use all observables in the fits and quote all errors at 1 a c.l., unless noted otherwise. Since 
our "fake" experimental data sets are perfect sets, the minimum of calculated equals zero 
and is thus not meaningful; only Ax^ calculated with respect to the best fit points has any 
physical meaning in the plots shown below. 
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lO^-' 10" lO^^ 10" 2 X 10" 5 X 10" lO^-' 2 x lO^^ 

niss [GeV] iriss [GeV] 

FIG. 6: Calculated distribution versus mss for 7 free parameters, P5 and mss = 10^"^ GeV (to 
the left) and mss = 5 x 10^^ GeV (to the right). 

Fig. (5) shows the allowed parameter space obtained in a MonteCarlo run for mo, M1/2, 
tan/3, Aq and mss for 7 free parameters, P5 and rngs = 5 x 10^"^ GeV. Shown are the allowed 
ranges of mo and M1/2 versus mss, as well as mo versus M1/2 and tan/3 versus Aq. On top 
of the 4 cMSSM parameters and mss in this calculation we allow the solar angle (^12) and 
the atmospheric angle (^23) to float freely within their allowed range. Errors on neutrino 
angles for this plot are taken from [63]. Plots for other points and/or different sets of free 
parameters look qualitatively very similar to the example shown in the figure. There is 
very little correlation among different parameters, contrary to the situation found in case of 
seesaw type-II and type-III [38]. Especially no correlations between mo, M1/2 and mss are 
found. However, there is some correlation between tan/3 and Aq, driven by the fact that 
m/jO alone can only fix a certain combination of these two parameters well. The correlation 
between tan /3 and Aq is slightly stronger than in the cMSSM case, due to the contribution 
of Aq in the running of slepton masses, see eq. (7). 

For our assumed set of measurements, mo and M1/2 are mainly determined by the highly 
accurate measurements of right slepton and gaugino masses of the ILC. Aq and tan /3 are fixed 
by a combination of the lightest Higgs mass and the lighter stau mass. LHC measurements 
help to break degeneracies in parameter space, but are much less important. We stress that 
the highly accurate determination of cMSSM parameters shown in fig. (5) is a prerequisite 
for determining reliable errors on m55'. ^ 

Fig. (6) shows calculated distributions versus mss for the same 7 free parameters as 
in fig. (5), P5 and mss = 10^'' GeV (to the left) and mss = 5 x 10^"* GeV (to the right). For 
the latter an upper (lower) limit of mss — 8 x 10^^ GeV {mss — 3 x 10^^ GeV) is found. For 
mss = 10^"^ GeV a clear upper limit is found, but for low values of mss the distribution 
fiattens out at Ax^ ~ 6.5. This different behaviour can be understood with the help of 
the results of the previous subsection, see fig. (2). For mss = 5 x 10^^ GeV, there exists a 
notable difference in some observables with respect to the cMSSM expectation, especially left 



^ We have checked this explicitly in a calculation using only LHC observables. 
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FIG. 7: Calculated allowed range of mss versus mss for 5 (left) and 7 (right) free parameters and 
P5. The two different error bars correspond to 1 and 3 a c.l. 



smuon and selection mass can no longer be adequately fitted by varying mo and M1/2 alone, 
without destroying the agreement with "data" for right sleptons and gauginos. Therefore 
both, a lower and an upper limit on mss exist for this point. The situation is different for 
mss = 10^^ GeV, for which the spectrum is much closer to cMSSM expectations. Larger 
values of mss are excluded, since they would require larger Yukawas, i.e. larger deviation 
from cMSSM than observed. Smaller values of mss, on the other hand, have ever smaller 
values of V^, i.e. come closer and closer to cMSSM expectations. For an input value of 
'^T'SS just below 77155' = 10^^ GeV there is then no longer any lower limit on 777.55, i-e. the 
data becomes perfectly consistent with a pure cMSSM calculation. In this case one can only 
"exclude" a certain range of the seesaw, say values of 77755 above a few lO^'^ GeV. 

One standard deviation is, of course, too little to claim an observation. We therefore 
show in fig. (7) A{mss) versus 77755 for 5 (left) and 7 (right) free parameters and P5 at 1 
and 3 a c.l. At 77755 = 10^^ formally a 1 sigma "evidence" could be reached, but at 3 cr c.l. 
the spectrum is perfectly consistent with a pure cMSSM. For larger values of 77755, however, 
several standard deviations can be reached. For the two largest values of 77755 calculated in 
this figure, a 5 a "discovery" is possible. 

Fig. (7) shows A(77755) for 5 and 7 free parameters. We have repeated this exercise for 
different sets of free parameters and 77755 = 5 x 10^^. Here, 5 free parameters correspond 
to the 4 cMSSM parameters plus 77755, 7 free parameters are the original 5 plus 9i2 and 
6*23. We have also tried other combinations such as 6 parameters: original 5 plus 6r and 8 
parameters, where we let all 3 neutrino angles float freely. Sets with larger numbers of free 
parameters are no longer sufficiently sampled in our MonteCarlo runs, so we do not give 
numbers for these, although in principle the calculation could allow also to let the neutrinos 
mass squared differences to float freely. Error bars are slightly larger for larger number of 
free parameters, as expected. However, since there is little or no correlation among the 
parameters, the differences are so small as to be completely irrelevant. 
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IV. SUMMARY AND DISCUSSION 



We have discussed the prospects for finding indirect hints for type-I seesaw in SUSY mass 
measurements. Since type-I seesaw adds only singlets to the SM particle content, only very 
few observables are affected and all changes in masses are small, even in the most favourable 
circumstances. Per-mille level accuracies will be needed, i.e. measurements at an ILC, before 
any quantitative attempt searching for type-I seesaw can hope for success, even assuming 
admittedly simplistic cMSSM boundary condtions. 

Our calculation confirms quantitatively that slepton mass measurements can contain in- 
formation about the type-I seesaw. Right sleptons are expected to be degenerate, while the 
left smuon and selectron show a potentially measurable splitting between their masses. If 
such a situation is indeed found, an estimate of m^^ might be derivable from ILC SUSY 
mass measurements. 

Above we have commented only on experimental errors. However, given the per-mille 
requirements on accuracy, stressed several times, also theoretical errors in the calculation of 
SUSY spectra are important. Various potential sources of errors come to mind. First of all, a 

1- loop calculation of SUSY masses is almost certainly not accurate enough for our purposes. 
We have tried to estimate the importance of higher loop orders, varying the renormalization 
scale in the numerical calculation. Changes of smuon and selectron mass found are of the 
order of the ILC error or even larger, depending on SUSY point and variation of scale. For 
the mass of the lightest Higgs boson it has been shown that even different calculations at 

2- loop still disagree at a level of few GeV [64]. Second, our calculation assumes a perfect 
knowledge of the GUT scale. Changes in the GUT scale do lead to sizeable changes in the 
calculated spectra for the same cMSSM parameters, which can be easily of the order of the 
required precision of the calculation and larger. In this sense, A(me;i) is an especially nice 
observable, since here the GUT scale uncertainty nearly cancels out in the calculation. In 
summary, if ILC accuracies on SUSY masses can indeed be reached experimentally, progress 
on the theoretical side will become necessary too. 

In our calculations, we have considered only SUSY masses. We have not taken into 
account data from lepton flavour violation, mainly because currently only upper limits are 
available. If in the future flnite values for /j — )■ ■ + 7 become available, it would be very 
interesting to see, how much could be learned about the type-I seesaw parameters in a 
combined flt. Including LFV one could maybe also allow for non-degenerate right-handed 
neutrinos in the flts. 

And, flnally, despite all the limitations of our study, we flnd it very encouraging that 
hints for type-I seesaw might be found in SUSY mass measurements at all. We stress again, 
that LFV and SUSY mass measurements test different portions of seesaw parameter space. 
For a more complete "reconstruction" of seesaw parameters, than what we have attempted 
here, both kinds of measurements would be needed. 
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